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Abstract 

Based on spin-density functional theory we calculate the electronic structure 
of a tunnel junction consisting of two magnetic Fe layers separated by an insulat- 
ing vacuum barrier selfconsistently. For the conductance the Landauer formula is 
evaluated in the ballistic limit as function of the magnetic configuration. Based 
on these conductances the tunnel magnetoresistance (TMR) ratio is obtained. We 
investigate the relation between TMR ratio and spin polarization of the electronic 
structure at the metal/insulator interface. 
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1 Introduction 



The discovery of Giant MagnetoResistance ^ in 1988 initiated a renewed interest in 
the phenomenon of spin dependent tunnehng. The first time Julhere in 1975 succeded 
to measure the spin polarized current in a tunnel junction composed of two ferromag- 
netic layers separated by an insulating barrier. Using improved experimental techniques 
it is now possible to produce high quality tunnel junctions showing a remarkable tunnel 
magnetoresistance 0, |] . A widely used model to interpret the results in terms of the po- 
larization of the ferromagnetic leads was proposed by Julhere [0. It was based on earlier 
studies of Meservey and Tedrow on tunneling between superconductors and ferromagnets 
@, 0, |]. The model assumed that the tunnel current is proportional to the product of 
the effective tunneling density of states of the considered spin direction for both ferro- 
magnetic leads. In this paper we want to focus on the ab initio description of TMR of a 
Fe/Vacuum/Fe junction using a supercell approach. We will discuss the properties of the 
local and total density of states and the influence of the finite thickness of the electrode 
on current polarization and TMR ratio. 



2 Transport 

To evaluate the tunnel current in an Fe/Vacuum/Fe system we use a supercell description. 
We consider unit cells with a fixed number of empty sites to simulate the vacuum barrier 
with a thickness of 0.55 nm (4 ML) and a variing number of Fe sites describing the 
electrodes. Based on density functional theory the electronic structure of the system 
was calculated selfconsistently using a Screened Korringa-Kohn-Rostoker (KKR) method 
B 0. This method is especially advantageous for prolonged supercells and we treated 
supercells with up to 108 atoms for this investigation. 

We consider ballistic transport perpendicular to the layers in the limit of zero bias. 
Starting from the Landauer formula the conductance through a barrier between two reser- 



voirs can be expressed by the following Fermi surface integral |TT], |T2|, [13 



G(n)=e2|5: 5{El-Ep)\u.^,ll (1) 

if we assume ideal transmission for all states at the Fermi energy Ep- A denotes the 
cross section of the junction, n • is the component of the Fermi velocity in current 
direction. The sum in Eq.(l) goes over all states k and over all spin directions cr, that 
means, we assume that the electron spin is conserved during the transmission process. 
The finger-print of magnetoelectronics is that the conductance chang GS clS cl function of 
the magnetic configuration from antiparrallel (AP) into parallel (P) under the influence 
of an external magnetic field. For this reason, we consider the electronic structure and 
the conductance of the superlattices in both magnetic configurations (see Fig. 1). The 
TMR ratio is finally defined by 

TMR = 1 - ■ (2) 
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Figure 1: Schematic picture of the considered structures, the directions of the Fe magne- 
tizations for P and AP configuration, the direction and the usual strength of the currents 
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Figure 2: Polarization of the local density of states at the Fermi level for a Fe(OOl) 
surface, Vc denotes the atomic sites above the surface 

3 Density of states and polarization 

Julliere assumed that the TMR ratio is proportional to the polarization of an effective 
density of states of the ferromagnetic leads. The question is which density is appropriate 
to satisfy this model. The polarization of any density of states is defined to be 



n 



(3) 



with rf being a spin- projected local density of states depending on the position in the 
unit cell or the total density of states for the whole system. Results for the polarization 
of the local densities of states of an Fe/Vacuum/Fe system with a vacuum barrier of 1.1 
nm (8 ML) are shown in Fig. 2. 

In Fig. 2 the local density of states at the Fermi level at the different sites in the 
supercell is shown. The Fe polarization is positive and bulk-like inside the Fe layers and 
changes sign at the Fe interface due to the existence of surface resonances. Inside the 
vacuum barrier the polarization remains negative and tends to zero for larger distances 
from the surface. 
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Figure 3: Dependencies on the Fe layer thickness of a) conductance for P and AP 
configuration, b) spin polarization Pq of P conductance G^, and c) TMR ratio; the 
horizontal lines indicate the average 



4 TMR 



Fig. 3a shows the corresponding conductance for P and AP configuration in dependence 
on the Fe layer thickness with a barrier consiting of 4 vacuum layers (0.55 nm). In 
the asymptotic region (for thick Fe layers) the conductance is a superposition of vari- 
ous oscillatory contributions around a constant value of 0.256 for P and 0.077 for AP 
configuration, respectively. All values of G are given in atomic units of the Rydberg sys- 
tem (Go = 1.752 ■ 10^'^ fl~^m~^). Characteristic oscillation periods are 2.3, 3.5 and about 
7.5 ML. These periods are related to stationary points on the bulk Fe Fermi surface and 
describe wave vectors at which an remarkable amount of states at Ep is quantized fl^, |15|. 

The Polarization P^ of the conductance in P configuration (Fig. 3b) behaves non- 
monotonic. The mean value in the asymptotic region is about -0.2. This is consistent with 
the negative polarization of the local density of states at the Fe surface and in the vacuum 
barrier region. MacLaren et al. |jl6[ obtained a positive polarization near 1 for a semi- 



infinite Fe/Vacuum/Fe junction. The difference should first be related to the structure of 
the barrier consisting of ASA potentials of constant height and secondly to the different 



4 



n(Ep) 



25 
20 
15 
10 
5 





• • P dn 

o o P up 
X X p 

o-o AP 



5 10 15 20 25 30 35 40 45 50 



#FeML 



Figure 4: Total DOS at the Fermi level of the Fe/Vc supercell in dependence on Fe layer 
thickness for P and AP configuration at a fixed barrier thickness of 4 ML (0.55nm); the 
horizontal fines indicate the Fe bulk values for both spin directions 

geometry of the leads. In our calculation the potentials are calculated self-consistently 
which leads to different shapes of the barrier potential. The barrier profile, however, can 
infiuence the transmission coefficients drastically as it is known from free electron models. 

The TMR ratio in dependence on the Fe layer thickness is shown in Fig. 3c. The 
horizontal line for thicknesses larger than 25 indicates the average of about 0.7. To 
elucidate the origin of the oscillations of conductance and TMR we investigated the total 
density of states of the periodic system at the Fermi level as a function of the Fe layer 
thickness (Fig. 4). Obviously, their is no difference for the density of states in P and AP 
configuration of the junction, that is, TMR is not caused by changes in the density of 
states. The enhancement of the density of states due to surface resonances dominates the 
minority band for thin Fe layers. For thicker Fe layers the spin-projected density of states 
tends towards the Fe bulk values indicated by the horizontal lines. The variations of the 
total density for large thicknesses arc very small and can not account for the oscillations 
of the conductances. That illustrates the fact that the conductance perpendicular to the 
layers is not caused by the huge number of states at the Fermi energy but by a few states 
with a large velocity in current direction. The relative position of these states with respect 
to the Fermi energy can change drastically with the Fe layer thickness due to quantum 
confinement and consequently causes the conductance and TMR oscillations. 

5 Summary 

The balfistic conductance was calculated for an Fe/Vacuum/Fe superlattice using the 
Landauer formula. Although we considered very large Fe layer thicknesses up to 50 
ML (about 7 nm) the conductance and the TMR ratio still oscillates. The origin of 
the oscillations is the finite size of the Fe electrodes. The calculated mean value of the 
conductance polarization Pq = —0.2 agrees with the polarization of the local density of 
states at the Fe surface. The calculated average of the TMR ratio is 0.7. 
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